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Using time-of-flight and triple-axis inelastic neutron spectroscopy, we determine spin wave exci- 
tations throughout the Brillouin zone for ferromagnetic manganites Lai-^Ca^MnOa (x = 0.25, 0.3) 
in their low temperature metallic states. While spin wave excitations in the long wavelength limit 
(spin stiffness D) have similar values for both compounds, the excitations near the Brillouin zone 
boundary of Lao.7Cao.3Mn03 are considerable softened in all symmetry directions compared to that 
of Lao.75Cao.25Mn03. A Heisenberg model with the nearest neighbor and the fourth neighbor ex- 
change interactions can describe the overall dispersion curves fairly well. We compare the data with 
various theoretical models describing the spin excitations of ferromagnetic manganites. 

PACS numbers: 75.30.Ds, 61.12.-q, 71.30.+h, 72.15.Gd 



I. INTRODUCTION 

The experimental investigation of spin dynam- 
ical properties in doped manganese perovskite 
Ax-xBxMxiOz, where A is the trivalent ion 
(La 3+ , Pr 3+ , Nd 3+ , etc) and B is the divalent ion 
(Ca 2+ or Sr 2+ ), is essential to the understanding of 
spin-spin interactions in these materials. At hole-doping 
level x « 0.30, these so-called colossal magnetoresistance 
(CMR) compounds exhibit an unusually large change 
in electrical resistance in response to a magnetic field 
and changes from a paramagnetic to a ferromagnetic 
state, 1 -! 2 -! 3 - The Mn 3d levels in the Mn 3+ /Mn 4+ mixed 
valent system, split by the oxygen octahedral crystal 
field to a lower energy t2 g triplet and a higher energy 
e g doublet, are filled according to the Hund's rule 
such that all spins are aligned on a given site by a 
large intra-atomic exchange coupling Jh- The basic 
microscopic mechanism responsible for the CMR effect is 
the double-exchange (DE) interaction, 4 - where ferromag- 
netism and electrical conductivity arise from hopping 
of the itinerant e g electrons with kinetic energy t from 
trivalent Mn 3+ to tetravalent Mn 4+ sites. In its simplest 
form, the Hamiltonian of a DE model can be described 
as a single band of itinerant e g electrons interacting 
with localized core spins by the Hund's rule exchange 
Jh- Since Jh is much larger than t, the kinetic energy 
of itinerant e g electrons is minimal when all electron 
spins are parallel, i.e., the ground state is a metallic 
ferromagnet. 

Although DE interaction is believed to be responsi- 
ble for the ferromagnetism and electron conductivity in 
CMR compounds, whether the magnetic excitations of 
such a model can be discussed in terms of an equivalent 
ferromagnetic Heisenberg model is still not clear. In the 
strong coupling limit (t/ 'Jh —> 0), Furukawa- has shown 
that the DE model can be mapped onto the Heisenberg 



Hamiltonian with only the nearest neighbor (NN) ex- 
change coupling. In this scenario, the magnitude of the 
exchange coupling J associated with ferromagnetic spin 
waves should scale with the Curie temperature T Cl kinetic 
energy t and doping x&£&& Experimentally, the initial 
measurements on Lao.7Pbo.3Mn03 suggest that a sim- 
ple NN Heisenberg model is sufficient to account for the 
entire spin wave dispersion relation 7 - and the exchange 
coupling obtained from such model also yields, to within 
15%, the correct T c of the compound. However, later 
experimental measurements indicate that spin wave ex- 
citations of most Ai-jB^MnOs manganites with x m 0.3 
renormalize near the zone boundary (ZB) with large soft- 
ening and damping ] 10 : 11 ! 12 : 13 ' 14 Furthermore, the spin 
wave stiffness constant D and the NN magnetic exchange 
coupling J are weakly dependent on T c and doping level 
x in metallic ferromagnetic (FM) manganitesi 15 ' 16 : 17 

It is now well established that a Heisenberg model with 
NN exchange coupling is insufficient to describe the dis- 
persion relation of CMR manganites, several possible mi- 
croscopic mechanisms have been proposed to address the 
unusual features of spin wave excitations near the ZB. 
First, realistic calculations based on the DE mechanism 
with consideration of the finite kinetic energy t and the 
effect of on-site Coulomb repulsion show that spin waves 
in the DE model do not map to the Heisenberg Hamil- 
tonian with simple NN exchange coupling^ This model, 
however, predicts a doping dependence on D which is 
not observed experimentally. Second, the observed ZB 
spin wave softening may be due to the conduction elec- 
tron band filling effect, where the existence of long-range 
magnetic interactions leads to the softening at the ZB. 19 
On the other hand, whether this approach is capable of 
explaining the observed spin wave broadening remains 
unclear. Third, large magnon-phonon interactions may 
give rise to the remarkable ZB softening along specific 
directions^ Fourth, the deviation of short wave length 



2 



magnons from the canonical Heisenberg form might orig- 
inate from the scattering of magnons by collective quan- 
tum orbital fluctuations, which could either be planar 
(x 2 — y 2 )-type orbitals associated with the A- type an- 



tifcrromagnetic (AF) orderin g u ' or rodlike (3z 
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orbital correlations related to C-type AF ordering. 22 De- 
pending on the actual orbital shape, the coupling be- 
tween charge and orbital-lattice will give rise to distinct 
doping dependence of the softening/broadening of the 
magnetic spectrum. Fifth, the randomness created by 
the substitution of the divalent ions for the trivalent ions 
in Ai- x B x Mn03 might be responsible for the anomalous 
spin wave softening ] 23 ' 24 Finally, the overlap between the 
magnon excitations and Stoner continuum in the metal- 
lic Ai- x B x Mn03 would cause softening and broadening 
of the magnon branch near ZBi 25 i 26 Although this single- 
band DE model with intermediate coupling can explain 
the softening/broadening in the low-T c compounds, re- 
markable similarities in systems with widely different T c 's 
indicate it is inadequate as the bandwidth of Stoner Con- 
tinuum is directly related to the T c 's. 

Given that there are so many possible models to ex- 
plain the ZB magnon softening, it is imperative to carry 
out systematic spin wave measurements and compare the 
results with predictions of various models. In a recent 
Letter—, we made such comparison for spin waves along 
the [£,0,0] direction and found that none of prevailing 
models can account for the data. In this article, we ex- 
pand our previous work and describe a systematic inves- 
tigation of spin wave excitations of the CMR manganitcs 
La .75Cao.25Mn0 3 (LCM025) and La .7oCao.3oMn0 3 
(LCMO30). Using reactor based and time-of- flight in- 
elastic neutron scattering (INS) techniques, we were able 
to map out the low temperature ferromagnetic spin wave 
excitations of LCM025 and LCMO30 throughout the 
Brillouin zone in all symmetry directions. In the long 
wavelength limit, spin wave stiffness of LCM025 and 
LCMO30 are 147 ± 3 and 169 ± 2 meVA 2 respectively, 
consistent with previous result a 16 ' 17 . At large wavevec- 
tors, we find that the dispersion relations of LCMO30 
are considerable more renormalized (softened) in all ma- 
jor symmetry directions compared to those of LCM025. 
The softening is well described by the introduction of the 
4 th NN ferromagnetic exchange coupling J4 [Fig. 1(b)], 
the ratio of J 4 /J x is about 19.5% in LCMO30 and 6.5% 
in LCM025. In section II, we describe experimental de- 
tails. Section III gives the data analysis and comparison 
with previous work. The conclusions are summarized in 
Section IV. 



II. EXPERIMENTAL DETAILS 



Bragg peaks. 27 LCM025 has a nominal hole doping level 
of x = 0.25, just above the metal-insulator transition 
concentration (x — 0.22). LCMO30 has a doping level 
close to optimal doping with highest T c . Our INS exper- 
iments were performed on the HET chopper spectrom- 
eter at the ISIS spallation neutron source, Rutherford- 
Appleton Laboratory, and on the HB1/HB3 triple- axis 
spectrometers at the High-Flux-Isotope Reactor (HFIR), 
Oak Ridge National Laboratory. The momentum trans- 
fer wavevectors q = (q x ,q y ,q z ) are in units of A -1 at 
positions (H,K,L) — (q x a/2ir, q y b/2ir, q z c/2ir) in recip- 
rocal lattice units (rlu), where a rs b rs c rs 3.87 A and 
3.86 A are the lattice parameters of the pseudocubic unit 
cells of LCM025 and LCMO30, respectively [Fig. 1(a)]. 
The samples were aligned in the (H, H, L) zone in both 
the ISIS and HFIR experiments. For the ISIS measure- 
ments, we use the HET direct-geometry chopper spec- 
trometer which has the 3 He filled detector tubes covering 
the scattering angles from 9-29° (PSD detectors covers 
from 2.5-7.9°). For the HFIR experiment, we use py- 
rolytic graphite as monochromator, analyzer, and filters, 
and the final neutron energy was fixed at Ef = 13.5 or 
14.7 meV. 
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FIG. 1: (Color online) (a) Crystal structure of 
Lai_ x Ca x Mn03 with distorted oxygen octahedra sur- 
rounding Mn ions, (b) Magnetic exchange interactions up 
to 4 th order between adjacent Mn ions. (c) Illustration 
of time-of- flight experiment. The energy- wavevector (E-q) 
region probed by the experiment intersects the spin wave 
dispersion surface, giving rise to spin wave rings projected on 
the scattering plane. For smaller energy transfer E, the ring 
centers at the zone center. For larger E, the ring centers at 
the ZB. 

For the HET measurements, the magnetic scattering 
intensities of the raw data were normalized to a vanadium 
standard. The scattering function 



We grew single crystals of LCM025 and LCMO30 us- 
ing the traveling solvent floating zone technique. The 
Curie temperatures of LCM025 (T c = 190 ± 1 K) and 
LCMO30 (T c = 238 ± 1 K) are determined from the elas- 
tic neutron diffraction on the (100) and (110) magnetic 



S(q,w) = 
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where fej and kt are the initial and final neutron wavevec- 
tors, respectively, the solid angle of scattering is fl, and 
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the energy transfer is Uu>. We used the program mslicbSS 
to visualize the q-u> data sets and to prepare the one- 
dimensional cuts along the high-symmetry spin wave di- 
rections, needed for further analysis using the program 
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TOBYFIT 



III. MODELING AND ANALYSIS 

The neutron scattering cross-section per formula unit 
(f.u.) for spin wave excitations is 



d 2 ndu 



7T. 9 2 /4 1 



( 2 ) 

where ("fro) 2 — 0.2906 barn, |F(q)| 2 is the magnetic 
form factor, g is the Lande factor (« 2), [n(u>) + 1] = 
1/[1 — cxp(— fiui/kBT)) is the detailed balance factor and 
\xb is the Bohr magneton. x"(q, w ) is the imaginary part 
of the generalized spin susceptibility which depends on 
the wavevector q and energy transfer Tuo. In the damped 
simple harmonic oscillator (DSHO) approximation, the 
normalized dynamical susceptibility x"(<l) w ) can be writ- 
ten as 

HI \ 4jUJUJ 

7T[(W 2 - LUq) Z + 4:{'JUJj Z \ 

where 7 characterizes the damping of the magnetic spins, 
ujq is directly associated with the spin wave dispersion 
relation. In the light damping limit, the intrinsic peak 
height A and width T of the spin wave excitation profiles 
are determined by the damping 7 and become 

A oc l/(27r 7 ),r oc 27. (4) 

The Hamiltonian for a Heisenberg ferromagnet is 

H = — — ^2 Jk&i ■ Si+fc, (5) 

i.k 

where Si denotes the magnetic moment at site i, and Jfc 
indicates the magnetic exchange coupling between neigh- 
boring sites. In an early study, a NN coupling J\ has been 
successfully employed to describe the entire spin wave 
dispersion relation E(q) for high-T c manganites. 7 Subse- 
quent measurements have shown the presence of ZB spin 
wave softening for all other manganites with dispersions 
being reproduced well by a Heisenberg Hamiltonian with 
higher order interactions . 10 ' 17 ! 22 Recently, we found that 
the introduction of the 4 th NN interactions J4 gives sat- 
isfactory description of ZB softening of the spin wave dis- 
persions for a wide range of doped manganites along the 
[£, 0, 0] direction.— We show here that such model also 
gives reasonable description of spin waves in all other 
symmetry directions. 



A. Results on Lao.7sCao.25Mn03 

LCM025 undergoes a ferromagnetic phase transition 
and becomes metallic below T c — 190 K. 27 To determine 



the dispersion of LCM025 at large momentum transfers, 
we measured its ferromagnetic spin waves with incident 
beam neutron energies (E{) of 32, 50, 75, 100, 125, 150, 
175, and 185 meV at T = 8.5 K (0.045 T c ) on HET. The 
sample was oriented such that either the [1,1,0] or [0,0,1] 
axis of the crystal is along the incident beam directions. 
For the neutron beam along the [1,1,0] axis of the crystal, 
we could get the dispersion relations along the [£, 0, 0] 
or [£,£, 0] directions. For the neutron beam along the 
[0,0,1] axis, the dispersion along the [£,£,£] direction can 
be obtained. 

Fig. [2] summarizes the spin wave excitations for inci- 
dent neutron energies of 50, 75, and 100 meV. Panels (a- 
c) show the two-dimensional color-coded contour plots of 
excitations in reciprocal space. Two rings of scattering 
are observed in these panels. The first and strongest of 
these two is centered at the (1,0,0) and corresponds to 
the intersection of the E-q region probed by the experi- 
ment, and the spin wave dispersion surface near the zone 
center. The second of these rings, centered at (1.5,0.5,0), 
corresponds to such intersections near the ZB because of 
a larger energy transfer E [Fig. 1(c)]. We cut the E-q 
data along the [£, 0, 0] direction, as shown in panels (d- 
f). The cut clearly shows two distinct peaks located at 
11.7 and 16.3 meV for E { = 50 meV [Fig. 2(d)]. These 
peaks gradually disperse outward with increasing inci- 
dent beam energy [Figs. 2(e), (f)]. As a function of in- 
creasing energy and approaching the ZB wavevector, the 
spin waves become broader in width and weaker in in- 
tensity [Fig. H(f)]. These results are consistent with our 
earlier measurements 

To obtain spin wave excitations along the [£,£,0] di- 
rection, we change the orthogonal viewing axes of the 
two-dimensional spin wave spectra and cut images along 
the [£, £, 0] direction. Fig. [3] summarizes the outcome of 
these cuts, it is clear that the spin wave peaks become 
broader and weaker near the ZB [Fig. |3fc),(f)]. 

We systematically cut the data along all high sym- 
metry directions for various incident beam energies. To 
obtain reliable exchange coupling constant J^'s and de- 
termine their directional dependence, we analyze data 
along one symmetry direction at a time with as many 
cuts as possible. The data are fit simultaneously using 
the dynamic susceptibility x"(w,q) described in Eqn. [3] 
with spin wave dispersion relation as 

foj<>(q) =A + 2S[J(0)- J(q)], (6) 
where A is the anisotropic gap and 

J(q)=^J fe e ic »-( Ri - R A (7) 

i 

Using the tobyfit non-linear least-square analysis 
program, we fit all the cuts by adjusting the peak ampli- 
tude A, damping term T and magnetic exchange cou- 
pling constant J^'s. The time-of- flight measurements 
do not provide information at small momentum trans- 
fers, but the data obtained using triple-axis spectroscopy 
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FIG. 2: (Color online) Spin wave excitations of LCM025 in 
the [£, 0, 0] direction. Panels (a-c) illustrate the 2D contour 
plots in the reciprocal space with incident energy of 50, 75, 
and 100 meV. The cross symbols indicate the FM zone center 
of (1,0,0). Panels (d-f) depict the corresponding scan profiles 
along the [£, 0, 0] direction shown by the arrows in the left 
panels. The typical cut width is 0.1 rlu. Solid lines are least 
square fits using the spin wave model described in the text. 
The energies associated with excitation peaks are labeled. 



showed negligible anisotropic spin gap near the zone 
center i 10 ! 11 ' 12 ! 13 ! 14 We therefore fixed A = in Eqn. [5] 
during the analysis. Following the results of our previous 
paper, we force Ji and J3 to be zero during the course 
of analysis, allowing only J\ and J4 to vary. 17 Table Q] 
shows the fitting results in LCM025 along three major 
symmetry directions of [£,0,0], [£,£,0] and [£,£,£]. We 
note that the value of exchange constants varies in a nar- 
row range for fitted results along different directions. For 
example, the NN exchange coupling Jj changes from 7.6 
to 8.4 meV and J4 varies from 0.26 to 0.58 meV. Using 
D = Vq w o(q)|q=o, where ojo(q) is the dispersion relation, 
we can calculate the spin wave stiffness constant D and 
found that D w 145 meV A 2 , consistent with the value 
obtained from low-g inelastic scattering measurement ji 7 - 
Since the fitting results along the three high symmetry di- 
rections are consistent, we included excitation data along 
all directions in the tobyfit. As shown in Table QJ the 
global analysis with as many as 38 data sets along all di- 
rections gives Ji = 7.83 ±0.06 and J4 = 0.51 ±0.03 meV, 
which yields the ratio of J 4 /J 1 = 0.065 ± 0.004 and stiff- 
ness constant D = 147.2 ± 2.6 meVA 2 . 

tobyfit gives only the fitting parameters J^'s, A and 
r for the spin wave dispersion relations. To actually con- 
struct a plot of E versus q of LCM025, we derive the 
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FIG. 3: (Color online) Spin wave excitations of LCM025 in 
the [£, £, 0] direction with incident energies of 50, 75, and 
100 meV. The viewing axis of panels (a-c) are rotated 45° 
with respect to those in Fig. [5] such that the scan profiles 
along [£, £, 0] can be obtained. The cross symbols indicate 
the FM zone center of (1,0,0). The dashed lines depict the 
cutting direction. 



energy at an individual wavevector E(q) by analyzing 
each single one-dimensional cuts as shown in Figs. d- 
f),^ . As shown in Fig. [4] the experimental data collapse 
nicely onto the solid line generated using Eqns. [JJ] and 
[7] with inclusion of the NN interaction Ji and 4 th NN 
interaction J4. Finally, we show that a fit to a purely 
NN exchange coupling J\ is not adequate. The disper- 
sion relations based only on J\ are plotted as the dashed 
lines, although these curves describe reasonable well the 
low-g data, they clearly deviate the data points at the 
ZB with higher energies. For example, the actual ZB 
energies along [£, 0, 0] and [£, £, 0] directions are 31.0 and 
62.5 meV, respectively. They are lowered by about 4 and 
8 meV from the prediction of the NN Heisenberg Hamil- 
tonian. It is noted that while the description using only 
the NN exchange coupling is not sufficient to characterize 
the entire dispersion relations in LCM025, this exchange 
constant J\ is similar to the value established in the INS 
study of Lao.7Pbo.3Mn03. 7 



B. Results on Lao.7Cao.3Mn03 

To determine the evolution of spin wave excitations 
of ferromagnetic CMR compounds, we also measured 
LCMO30 using HET. For these measurements, we used 
incident beam energies of 50, 75, 100, 125, 150 and 185 
meV. Fig. [5] and Fig. M show raw data for LCMO30 
with the same incident energies as those of LCM025 
in Figures and [3] It is clear that LCMO30 has a 
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Direction Ji(meV) J4(meV) J4/J1 D(meVA 2 ) No. of data sets \ 2 

[£,0,0] 7.56(7) 0.58(4) 0.076(6) 147.1 ±3.5 11 1.45 

[f,£,0] 7.99(9) 0.44(4) 0.055(6) 145.6 ± 3.9 13 1.81 

[£,£,£] 8.39(14) 0.26(8) 0.032(10) 140.7 ±6.6 14 2.09 

[£,0,0] + [£,£,0]± [£,£,£] 7.83(6) 0.51(3) 0.065(4) 147.2 ± 2.6 38 1.90 

[g, 0, 0] + [g, £, 0] + [j, g, g] 8.75(3) O(Fixed) 130.3 ± 0.5 38 2.32 



TABLE I: Fitting parameters of the exchange coupling constant J\, J4, the ratio of J4/J1 and spin wave stiffness constant D in 
LCM025. Each data set is one dimensional scan profile as shown in Fig.[2jd-f) and consists of one or more magnetic excitation 
peaks. The results with only the NN interaction Ji are listed in the last row for comparison. 




FIG. 4: (Color online) The spin wave dispersion curves for LCM025. The solid points are data collected at ISIS and the open 
points at low energies are collected at HFIR. The dash (red) curves are the fits using only the NN interaction Ji. The solid 
(blue) curves are the fits combining the NN interaction Ji and the 4 th NN interaction J4 (see text). 



softened dispersion near the ZB. As shown in Fig. [H^f), 
the peak height at q = [1.4,0,0] is only half of that at 
q = [0.7,0,0], compared with 80% for that in LCM025 
[Fig. 2(f)]. From Eqn. 4, the decrease in peak height 
at a similar wavevector indicates a larger damping term 
7. In addition, the spin wave ring away from (1.5,0.5,0) 
[Fig. [5(c), (f)] collapses much faster than in LCM025, 
as the incident beam energy increases to 100 meV. Ta- 
ble [II] lists the quantitative fitting parameters using TO- 
BYFIT. Much reduced J\ and considerable enhanced J4 
are observed along different high symmetry directions for 
LCMO30. Similar to LCM025, the softening near the 
boundary shows little symmetry directional dependence. 
For the global fit, we obtained J\ — 6.36 ± 0.03 and 
J 4 = 1.24 ± 0.02 meV. The ratio of J 4 /J x reaches 20%, 
a value much larger than that in LCM025. The stiffness 
constant is calculated to be D = 169 ± 2 meVA 2 , agrees 
well with the early triple-axis scattering result . 16 > 31 

The dispersion curves covering the entire Brillouin 
Zone of LCMO30 are shown in Fig. [7] with the solid lines 
representing the fits using J\ and J4, while the dashed 
lines representing the fits using only J\. Inspection of 
Fig. [7] reveals that, as in LCM025, the inclusion of the 
J4 gives a better fit to the data. The normalized chi- 
squared, x 2 > including J4 is a factor of 2 less than that 



without it. The renormalization of spin wave disper- 
sions at large wavevectors is better seen in the [£,0,0] 
and [£,£,0] directions, where the boundary energies are 
lowered to 25 and 52 meV, respectively. They are nearly 
9 and 15 meV lower than the expectation from a simple 
NN Heisenberg Hamiltonian. We note that fits at larger 
energies (E > 70 meV) are not satisfactory, possibly be- 
cause uncertainties associated with the determination of 
a already damped magnetic excitation at high energies. 

Having determined the dispersion relations for 
LCM025 and LCMO30 using the model described in 
Eqn. [5] with inclusion of J\ and J4, we now consider the 
other important aspect of the spin dynamics, the intrin- 
sic width r(q). This lincwidth T(q) is calculated from 
Eqn. 4, where 7 is obtained from the best fit to Eqn. [3] 
It is directly associated with the relevant damping mech- 
anisms and reflects how the quantized magnons interact 
with other scattering processes. The wavevector depen- 
dence of r(q) along different directions is shown in Fig. [5] 
Note that there is marked difference between LCM025 
and LCMO30; the widths of the latter are always larger, 
indicating more damped excitations in LCMO30. In ad- 
dition, the momentum evolution of T(q) does not show 
any anomaly across the Brillouin zone, appear to be iso- 
topic and increase drastically near the ZB. For example, 
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Direction 


Ji(meV) J 4 (meV) J4/J1 


D(meVA 2 ) No. of data sets 


x 2 


K,o,o] 


6.01(5) 1.29(4) 0.215(8) 


166.6 ±2.8 


13 


1.71 




6.52(5) 1.18(4) 0.181(6) 


167.4 ± 2.5 


22 


1.91 




6.56(9) 1.18(6) 0.179(11) 


167.8 ±4.8 


15 


1.54 


[£,0,0]±[£,£,0]± 6.36(3) 1.24(2) 0.195(4) 


168.6 ± 1.7 


50 


1.74 


K,o,o] + K,€,o] + K,f,€] 


8.43(3) O(Fixed) 


125.7 ±0.5 


50 


3.25 



TABLE II: Fitting parameters of exchange coupling constant Ji, J4, the ratio of J4/J1 and spin wave stiffness constant D in 
LCMO30. See Table U for additional information. 
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FIG. 5: (Color online) Spin wave excitations of LCMO30 in 
the [H,0,0] direction. Panels (a-c) illustrate 2D contour plots 
with incident energies of 50, 75, and 100 meV. Solid lines are 
least square fits. See Fig. [2] for additional information. 



r(q) reaches around 20 meV for E greater than 60 meV 
along the direction. This is unexpected from a 

classic, cubic Heisenberg ferromagnet, where the mag- 
netic excitations near the ZB are well resolved^ The in- 
trinsic widths in both LCM025 and LCMO30 are much 
broader, indicating that one or more decay mechanisms 
play an important role in these CMR compounds. 

We are now in a position to compare our data to var- 
ious possible mechanisms of the ZB softening. First, we 
want to comment on the role of disorder or randomness 
which is naturally present in these compounds. Motome 
and Furukawa3£ have pointed out that disorder in CMR 
materials will cause anomalous broadening and/or anti- 
crossing in the spin excitation spectra. In this scenario, 
the one-electron bandwidth is proportional to the average 
ionic-size at La/Ca site (f = J^. a^r, where Xi is the frac- 
tional occupancies of A-site species, is the individual 
radius). In a previous paper, we have characterized the 
softening of dispersion relation in a series of doped man- 
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FIG. 6: (Color online) Spin wave excitations of LCMO30 in 
the [H,H,0] direction with incident energies of 50, 75, and 
100 meV. See Fig. [2] for additional information. 



ganites near x = 0.30 along [£, 0, 0] direction^ We found 
the renormalization near ZB has little dependence on the 
average ionic size at A-site. For LCM025 and LCMO30, 
f become 1.207 and 1.205 A, respectively, showing very 
little variation of oxygen octahedron distortion surround- 
ing the Mn-ions. On the other hand, the mismatch be- 
tween La and Ca ions will cause a quenched disorder in 
the system, which can be qualitatively characterized by 



^ 2 = £*( 



Xir; 



2) ,33,34 Base on thiSj CT 2 jg 2.43 x 10 " 



for LCM025 and 2.72 x 10~ 4 for LCMO30. Clearly, the 
differences in quenched disorder between LCM025 and 
LCMO30 are rather small and cannot account for the 
dramatic change of spin wave spectra near the ZB. We 
therefore conclude that the disorder effect does not play 
an important role in this doping range of LCMO. 

Second, the effect of magnon-phonon coupling may be 
the microscopic origin of the observed magnon broad- 
ening and damping. Dai and co-workers suggest that 
the interaction between optical phonon and spin wave 
branches may lead to the broadening of spin wave spec- 
trum in a number of doped manganites near x — 0.30. 12 
In this picture, a dispersionless optical phonon branch 
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FIG. 8: The wavevector dependence of intrinsic excitation 
widths r(q) along three high symmetry directions in LCM025 



and LCMO30. The g-dependence of T(q) along [£, 0, 
LCMO30 obtained from triple-axis measurement (Ref. 
is displayed for comparison. 



for 
II) 



with energy around 20 meV goes across the whole Bril- 
louin zone and interacts with the magnon branch, causing 
softening and a substantial increase of the magnetic ex- 
citation linewidths. The measurements of Dai et al. were 
carried out with unpolarized neutrons, and it was diffi- 
cult to separate the magnetic scattering from the purely 



lattice excitations. Recent measurements by Fernandez- 
Baca et al£^ using polarized neutron scattering tech- 
niques confirmed that the ZB spin waves of LCMO30 in 
the [£,0,0] and [£,£,0] directions are considerable broad 
and have energies lower than those expected from the 
approximation to the NN Heisenberg Hamiltonian, al- 
though the softening in [£, £, 0] is less as severe than orig- 
inally reported?^ While the magnon-phonon interaction 
mechanism seems to explain the broadening of the spin 
waves, it is not clear if it would fully account for the 
magnitude of the observed softening in LCMO30 at the 
zone boundary. Furthermore, this mechanism may not 
be relevant to the general case of the manganites as the 
observed spin waves in the Smo.55Srn.45Mn03 are soft- 
ened to around 15 meV (which is below the 20 meV op- 
tical phonon) at the ZB along [£, 0, 0] direction, with no 
evidence of broadening^ 

Finally, the distinct feature of ZB softening might 
be a consequence of the e s -band filling in the half- 
metallic region, 12 Solovyev and Terakura suggested that 
the canonical double exchange is no longer appropriate 
as soon as holes are doped into the system. Longer range 
FM interactions lead to the softening at the zone bound- 
ary and contribute to the increase of the stiffness constant 
Z?, 19 However, the details of the realistic electronic struc- 
ture are important and may significantly modify the anal- 
ysis, particularly if the effects of ti g electrons are taken 
into account. For example, the x dependence of mag- 
netic interactions is substantial modified by the change of 
DE interaction contributed by ti g electrons and D might 
even decrease with x. Such complex scenario prevents 
us to make meaningful comparison with the D's yield 
experimentally. It is also unclear whether this model 
can explain the commonly observed damping/broadening 
of magnetic excitation near the ZB. We would like to 
point out that, in the tight-binding approximation used 
by these authors, only the exchange couplings along the 
Mn-O-Mn chain (Ji, J4, J§, and J15 as defined in Ref.[l9|) 
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would bring appreciable contributions to the spin wave 
dynamics. Our work, and that of Endoh, — show that Ji 
and J4, are the only exchange constants contributing to 
the ZB softening. There is close correlation between Jj's 
and the orbital polarized states* 2 ^ J2 would be enhanced 
by (x 2 — y 2 )-type orbitals and J4 is enhanced of (3z 2 — r 2 )- 
type orbitals. It is surprising that a large increase of 
Ji/ Ji occurs despite of the small change in nominal hole 
doping. This might indicate a drastic modification of 
(preformed) orbital correlations which favor the overlap 
between neighboring Mn ionsi^l We hope that the results 
presented here will help to stimulate further experimen- 
tal and theoretical investigations leading to a complete 
understanding of the magnetic dynamics in those CMR 
materials. 



IV. CONCLUSIONS 

In summary, we have performed a systematic study 
on the spin wave dynamics in the CMR manganites 
LCM025 and LCMO30 using inelastic neutron scattering 



techniques. We find that both systems display consider- 
able spin wave renormalization along all high symmetry 
directions. Consistent with early measurements, which 
entirely focus on the [1,0,0] direction, we find that the 
dispersion relations can be phenomenologically analyzed 
using the NN interaction J\ and the 4 th NN interaction 
Ji. The introduction of J4 lowers J\ and therefore low- 
ers the ZB energy. As a result of this the systems with 
a larger J4/J1 ratio exhibit a larger SW softening at the 
ZB. The possible mechanisms responsible for such bound- 
ary softening are also discussed. 
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